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The two-dimensional layered perovskite Sr2lr04 was proposed to be a spin-orbit Mott insulator, 
where the effect of Hubbard interaction is amplified on a narrow Jeff = 1/2 band due to strong spin- 
orbit coupling. On the other hand, the three-dimensional orthorhombic perovskite (Pbnm) SrIrOs 
remains metallic. To understand the physical origin of the metallic state and possible transitions 
to insulating phases, we construct a tight-binding model for SrlrOs- The band structure possesses 
a line node made of Jeff = 1/2 bands below the Fermi level. As a consequence, instability towards 
magnetic ordering is suppressed and the system remains metallic. This line node, originating from 
the underlying crystal structure, turns into a pair of three-dimensional nodal points on the introduc- 
tion of a staggered potential or spin-orbit coupling strength between alternating layers. Increasing 
this potential beyond a critical strength induces a transition to a strong topological insulator, fol- 
lowed by another transition to a normal band insulator. We propose that materials constructed with 
alternating Ir- and Rh-oxide layers along the (001) direction, such as Sr2lrRh06, are candidates for 
a strong topological insulator. 

PACS numbers: 73.20.-r, 71.30.+h, 71.70.-d 



I. INTRODUCTION 

Recently, there has been intense theoretical and exper- 
imental studies on iridium oxides, dubbed iridates, due to 
the possibility of realizing spin liquid^ ^, topological band 
insulatoi^ ^, and topological Mott insulatoi^^H^. Unlike 
3d transition metal oxides, 5d materials are expected to 
be less correlated because of extended orbitals. However, 
due to strong spin-orbit coupling (SOC), the Jeff = 1/2 
band formed out of t2^-orbitals is rather narrow, lead- 
ing to stronger correlation effect of the on-site repulsive 
Hubbard interaction U. In particular, pyrochlore iridates 
have been widely studied due to their geometrical frustra- 
tion, and novel topological phases have been theoretically 
proposed.™lini 

The perovskite iridates, on the other hand, have at- 
tracted less theoretical attention because of the lack of 
frustration. Due to this non-frustrated lattice structure, 
an antiferromagnetic insulator is expected when the Hub- 
bard U is large compared to the kinetic energy scale 
as measured by the bandwidth. It was reported that 
the SOC and Hubbard interaction are about the same 
magnitude in two-dimensional (2D) layered perovskite 
Sr2lr04^^ and a combination of the SOC and Hub- 
bard interaction results in a canted antiferromagnetic in- 
sulator (AFI). This phenomena was well described by a 
strong Hubbard U limit in Ref. TT, While Sr2lr04 has 
a canted AFI ground state, three-dimensional (3D) or- 
thorhombic perovskite SrIrOa remains metallic.'^^H^ 

In this paper, we investigate the band structure of or- 
thorhombic perovskite with Pbnm symmetry applicable 
to SrIrOa using a tight binding approach. We show that 
there is a line node near the Fermi level originating from 
the underlying lattice structure. Thus, SrIrOa is metal- 
lic with small Fermi pockets and becomes a perfect 3D 
semi-metal with nodal dispersion when the SOC becomes 



stronger. When a staggered potential or SOC between 
alternating layers is introduced, the nodal-line is gapped 
out leaving a pair of 3D nodal points. The system under- 
goes a transition to a strong topological insulator beyond 
a critical strength of the potential. A further increase of 
the potential turns the system into a trivial band insula- 
tor. The phase diagram is illustrated in Fig. [l] and the 
nodal band structure along with the mechanism for the 
transition are presented below. 



II. SYMMETRY OF LATTICE AND BAND 
STRUCTURE 

The Ruddlesden-Popper series for the strontium iri- 
dates Sr^+iIrnOsn+i is a family of materials with Ir 
atoms surrounded by an octahedral environment of oxy- 
gen. The octahedral crystal field splits the d electrons to 
low-energy t2g and high energy e^, levels. The large spin- 
orbit coupling further splits the t2^ electrons into low- 
lying Jeff = 3/2 states and higher-lying Jeff = 1/2 states. 
In these compounds, the oxidation state of the Ir atoms 
is Ir^+ leaving five d-electrons on the outer shell. When 
the SOC is larger than Hund's coupling, the Jeff = 3/2 
band is completely filled while the Jeff = 1/2 band is 
half-filled, leading to an expected metallic state based on 
band theory. 

Among the quasi-2D perovskites, Sr2lr04 is a single 
layered material where the iridium-oxide layers are sep- 
arated by strontium atoms^^ and the octahedral envi- 
ronment of oxygen around the iridium atoms is rotated 
around the z-axis in a staggered pattern. Sr3lr207 is 
a bilayer version of the single layered materiaPSl^ with 
opposite rotation of the octahedra in the second layer. 
Both exhibit a canted antiferromagnetic insulating phase 
due to the effect of the Hubbard interaction U On 
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FIG. 1. [Color online] Phase Diagram as a function of the 
staggered potential between layers, m. The product of parity 
eigenvalues for the filled bands are denoted by (=b) for each 
time reversal invariant momentum (TRIM) point. When m is 
finite, the line node becomes a nodal point located between U 
and R. This node gets shifted toward R (as indicated by the 
red arrow) as m increases until it reaches md, after which a 
full gap appears. Beyond md, the system turns into a strong 
topological insulator with indices uq; (z^iz^2^3) = 1; (001). At 
mc2, the gap closes at the Z point and the system becomes a 
topologically trivial band insulator for m > mc2- 



the othe r han d, SrIrOs has an orthorhombic perovskite 
stmctur j^^ l -^^ l where the octahedra are rotated around the 
(110) direction, in addition to the rotation about the z- 
axis. This yields a 4-atom unit cell with the space group 
of Pbnm shown in Fig. |2j We model the different oc- 
tahedral environment as a rotation around the z-axis by 
an angle of followed by another rotation around the 
new (llO)-direction by an angle of The colors rep- 
resent Ir atoms with different oxygen environment and 
corresponds to Blue = {0^(p)^ Red = (—6>, —(/:?), Yellow 
= (6>, —(p)^ and Green = (— 6>, (p). 



A. Symmetries of Pbnm 

The Pbnm structure, in addition to inversion, has the 
following symmetries: a b-glide, an n-glide and an m re- 
flection. Here, the a-axis is given by the (llO)-direction, 
the b-axis by the (llO)-direction and the c-axis is con- 
cordant with the z-axis. 

The b-glide is represented by a reflection of the a-axis 
(mirror b-c plane at a = 1/4) followed by a translation 
along the b-axis. The n-glide is represented by a reflec- 
tion of the b-axis (mirror a-c plane at b = 1/4) followed 
by a translation along the a-axis and another one along 
the c-axis. The m-reflection is represented by a reflection 
of the c-axis (mirror a-b plane at c = 1/4). These symme- 
try operators are denoted by XI^, 11^ and 11^ respectively 
The effect of these operations on the Ir atom position, the 
crystal momentum, the orbitals and the spin is given in 
Table |lj Here the r and u Pauli matrices represent the 
two sublattice indices, with r for the in-plane sublattice 
(i.e. blue and red or yellow and green) and u for the 
sublayer one (i.e. blue and yellow or red and green). 




FIG. 2. [Color online] Crystal structure for SrIrOa from a 
side view. The different colors of blue (B), red (R), yellow 
(Y) and green (G) represent Ir atoms with different oxygen 
environment. These 4 Ir atoms form a unit cell 



TABLE I. Effect of Pbnm symmetry operations on states. 
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Since the J = 1/2 basis can be written as 

I = ±1/2) = -^{\dyz^ Ts) ± i\d:,^,Ts) ± \d^y, ±s)), 

the orbital and spin (s) transformation can be written 
for the pseudospin J = 1/2 using the a Pauli matrices. 
Doing so, we can summarize the Pbnm symmetry oper- 
ations as 

i 



H77}, = io'zl^xkm') 



(1) 



where /c^, kn and km are the operators that act on the 
crystal momentum as shown in the third column of Table 

m 



B. First principles calculations 



Density functional theor^E^l^ calculations including 
Hubbard U interaction and SOC have been performed 
using the full-potential linearized augmented-plane-wave 
(FP-LAPW) method as implemented in the elk code.l^ 
We used the Ceperley- Alder form of the local-density ap- 
proximation (LDA)^^ as parametrized by Perdew and 
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Zungei^^ and used the "around mean field" (AMF) 
scheme^^ for the double-counting correction to the 
electron-electron interaction. To see the effects of reduc- 
tion or enhancement of the SOC, we simply multiply the 
corresponding SOC term in the Hamiltonian with a factor 
a; a = means no SOC and a = 1 means atomic SOC as 
implemented by the elk code. In the FP-LAPW method, 
real space is divided into spheres around the atoms (muf- 
fin tins) and interstitials elsewhere. In our calculations, 
the mufhn-tin radii determined automatically by the elk 
code are 1.86, 2.08, and 1.51 a.u. for strontium (Sr), irid- 
ium (Ir) and oxygen (O), respectively. The basis set we 
used consists of APW functions with angular momentum 
/ up to eight and plane waves with cut-off energy equal 
to 231.3 eV. The number of empty states in the basis 
set in the second variational step was 10. Brillouin zone 
integrations were performed using a 3 x 3 x 3 grid, which 
is equivalent to using 10 points in the irreducible part of 
the Brillouin zone. We only used U for the 5d orbit als of 
the Iridium atoms. 

The computed band structures for U = 2.0 eV with 
a = and 2 are shown in Fig. |3] Note that SOC 
significantly changes the band structure. In particular 
near the Fermi energy, the bands are mostly composed of 
Jeff = 1/2 states for a = 2. The band structure near the 
Fermi energy remains qualitatively similar from U = to 
U = 2 eV, and therefore only the U = 2 eV case is shown 
here. An important characteristic is the presence of the 
node near the U point revealing the semi-metallic prop- 
erty of SrIrOa. To understand this feature and possible 
transition to topological insulating phases, we derive a 
tight-binding model, described in the next section. 
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FIG. 3. LDA calculation with Hubbard [/ = 2 eV and SOC 
tuning parameter a — 0(a) and a — 2.0(b). The band struc- 
ture remains qualitatively similar from [/ = up to [/ = 2 eV 
while SOC has a significant effect. See the main text for the 
comparison to tight binding model. 



III. TIGHT BINDING MODEL 

A tight-binding model is constructed, taking into ac- 
count only the nearest and next-nearest neighbor Ir-Ir 
direct hoppings between Jeff = 1/2 orbitals, using the ba- 
sis = (CkBt5Cki?t'^kyt'^kGt5CkBi,Cki?|,Cky|,CkGi) 

where {B,R,Y,G} = {Blue, Red, Yellow, Green} de- 
scribed above, and (t^i) represent J^^ = 1/2 and —1/2, 
respectively. Note that first principles calculations show 
that all the Jeff = 3/2 states are sufficiently below the 
Fermi energy and the Jeff = 1/2 states occupy bands 
near the Fermi energy. Thus the current tight binding 
model of Jeff = 1/2 orbitals is expected to provide rele- 
vant bands near the Fermi level. 

In such a Jeff = 1/2 basis, the single-particle Hamilto- 
nian is written as 

= Ak + Re{El)r, + Im{El)a,ry + {E^)u, 

^Re{E^)cryVzTy + Im{E^)crxh'zTy 
-^Re{E^'')GyyyT^ + Im{E^'')G:^UyT^ 
-^Re{E^'')ayUyr^ + Im{E^'')a^Uyr^, (2) 

where the superscripts p, and d represent hopping be- 



tween Ir atoms within a layer (e.g., B to i?), between 
layers (e.g., B toY), and between next-nearest neighbors 
(e.g, B to B within a layer or 5 to G between layers), 
respectively, and o denotes hoppings between J^^ = 1/2 
and —1/2.'^ The coefficients are given by 

Ak = A + txy cos(kx) cos(ky). 

El = {2t - it'){cos{kx) + cos(ky)), 

E^ = 2t cos(kz), E^"" = (1 - iK cos(kz), 

E^ = tc^(cos(kx) + cos(ky)) cos(kz) 

+zt[^(sin(kx) + sin(ky)) sin(kz), 
(^?pCOs(kx)+t^^cos(ky)) 
-i {t^p cos(kx) + cos(ky)) , 
= t2(sin(ky) + i sin(kx)) sin(kz). (3) 

where A denotes the SOC, t nearest neighbor intra-orbital 
hoppings, txy next-nearest neighbor in-plane hopping for 
xy orbitals, and t' the hopping between ID-orbitals (xz 
and yz) in the plane, td and t'^ are originated from intra- 
orbital diagonal hopping of next-nearest neighbor atoms 
from different layers (e.g., B to G). comes from ID 
to xy hoppings between two adjacent atoms on different 
planes whereas tf^ and represent nearest neighbor in- 
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plane hoppings between ID and xy orbital. Finally, 
arises from out-of-plane next-nearest neighbor ID to xy 
orbitals hoppings. It is important to note that t' ^ t'^ , t^, 
tfp, and tJJ vanish when there is no rotation and/or 
tilting of the oxygen octahedra. 

Note that the symmetry operators of the Pbnm struc- 
ture described in section III Al all commute with the 
Hamiltonian given by Eq. |2j 

The hopping parameters are obtained using the Slater- 
Kosteil^ method and determined by {ta^t^^ts) for the 
nearest neighbor hoppings and {t'^^t'^^t'^) = a^{tcj,t^,ts) 
for the next-nearest neighbor hoppings. Thus, the only 
free parameters are {tcr^t-j^^ts) and a, and the unit of 
energy is set as l^^^l = 1 throughout the paper. All 
the other parameters are determined by the rotation 
and tilting angles of the Pbnm structure. We choose 
{tcr^tj^^ts) = (3/2,-1,1/4) and a = 0.4 for the values 
of the free parameters to match the LDA band struc- 
ture shown in Fig. [sj^b). With the angles {0,(f) = 
(11^12°)l^ this yields (t, t^ t,^, t?, t^, t'^, tf^, t^^, t^) = 
(-0.6,-0.15,0.13,-0.3,0.03,0.1,0.3,0.06) for the tight- 
binding parameters. 

The band structure is shown in Fig. 4(a) While the 
detailed dispersion is sensitive to the rotation and tilting 
angles, and to the choice of a parameter set, the robust 
feature is a line node near the U point in the Brillouin 
zone, consistent with first-principles calculations. This 
line node, which is underneath the Fermi level as shown in 
Fig. 4(a)[ originates from the reflection symmetry of the 
crystal structure at the z = ^ and ^ = | planes present 
in the Pbnm space group (Ilm)- Without t^, there is 
a node that is 8- fold degenerate at the U = (f , — f , f ) 
point. With a finite t^, a gap opens at the U point, and 
the nodes are pushed away from the U point in the XURS 
plane. It was found that an elliptical line node, shown 
as the red line in Fig. |4(b)[ appears in the XURS plane. 
The line node obtained by expanding the dispersion near 
the U point up to linear terms is given by 
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FIG. 4. [Color online] (a) Underlying band structure of or- 
thorhombic perovskite SrlrOs obtained by a tight binding 
model including only Ir-Ir direct hoppings. Note that the 
band structure here is qualitatively similar the the LDA cal- 
culation of Fig. [s] especially near the Fermi surface, (b) 
Brillouin zone with all TRIM points labeled along with the 
product of the parity eigenvalues of the filled bands. The red 
line in the XURS plane corresponds to the line node in the 
system. 



(4) 



where, p+ and are defined such that = \ +p+ 

= — I + p+ — p_ and = f +^2- The plane of the 
ellipse appears at p_ =0.^^ 

Due to this robust feature of nodes originated from 
the combination of SOC and Pbnm symmetry of the lat- 
tice, SrIrOa has small Fermi pockets with small density 
of states near the Fermi energy. As a consequence, it is 
expected that instability toward any magnetic ordering 
is suppressed. Indeed a larger value of the local Hubbard 
interaction is required for a canted antiferromagnetic or- 
dering, when compared to the quasi-2D perovskite iri- 
dates. As shown in Fig. [3j the system remains metallic 
for U up to 2 eV. Since the line node is protected by lat- 
tice symmetries and is not affected by moderate values of 
Hubbard interaction, perturbations that break the sub- 
layer reflection symmetry (Hrn) at 2: = | and | planes 
are required to open a gap along the U to R line and turn 



the system into an insulating phase. Below, we discuss 
one such perturbation and analyze the resulting phases. 



IV. TOPOLOGICAL INSULATOR 

The mirror symmetry of the Pbnm space group (Hm) 
commutes with every term of the Hamiltonian. In order 
to break such symmetry, we introduce a staggered sub- 
layer potential. Such a staggered chemical potential is 
added to the Hamiltonian in Eq. [2] as a mass term of 
the form mv^^ which anticommutes with H^, and hence 
breaks this symmetry. This term can also arise when the 
SOC alternates between the layers (B-R vs Y-G planes 
along 2;-axis as shown in Fig. 2). As the mass term is 
introduced, the degeneracy of the line node is lifted ex- 
cept for a pair of points along the U to R line which 
gets shifted towards the R point as shown by the red ar- 
row in Fig. [1] Along this line, only two terms of the 
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Hamiltonian are non-zero, only one of which commutes 
with muz^ having the effect of pushing the node instead 
of hfting it completely. At m = md, the node reaches 
the R point, and the system changes into an insulating 
phase for m > rrici as shown in Fig. [T] The insulating 
phase, above this critical value of the mass term, is classi- 
fied as a strong topological insulator characterized by the 
Z2 topological indices. ^^ ^^ The product of parity eigen- 
values of the filled bands at each TRIM points is shown 
in Fig. [1] and the topological indices for this phase are 
z^o; (^i^2^3) = 1; (001). Notice that the product of the 
parity eigenvalues at the R point has flipped sign at md . 

Further increase of the mass term leads to another 
transition at mc2 where the bands invert at the Z point, 
and the product of parity eigenvalues switches sign ac- 
cordingly as shown in Fig. [l] Beyond mc2, the system 
becomes a trivial band insulator. For the set of parame- 
ters we have chosen, the values of md and mc2 are found 
to be approximately md = 0.27 and mc2 = 2.4, forming 
a large window of strong topological insulating phase. 




FIG. 5. [Color online] The edge state calculation with a finite 
slab of material having periodic boundary conditions along 
the X and y directions and open boundary condition along 
the z direction. From the Fu-Kane criteria, we expect to see 
an odd number of crossing of the edge state between F and 
any of the other three TRIM points in the surface Brillouin 
zone. Here, we compute the band structure of such a system 
between F — (0, 0) and S — (tt, 0) for m — 0.5. 

To confirm the above Z2 analysis, we compute the sur- 
face state of a fini te sla b to verify the presence of pro- 
tected edge stdites^^^ The band structure calculation 
using periodic boundary condition along the x and y axis 
and open boundary conditions along the z axis is shown 



in Fig. [5] The red lines represents the surface states 
and the calculation was done using a Blue-Red layer as 
one edge, and a Yellow- Green layer as the other edge. 
Because of the sign difference of the mass term on these 
two layers, the surface states coming from these two edges 
have opposite velocity. 

V. DISCUSSION AND SUMMARY 

Based on the results of band structure obtained for 
SrIrOa with Pbnm space group symmetry, a staggered 
potential in alternating layers can turn the semi-metal 
state into a strong topological insulator. Such a potential 
can be obtained if the Ir atoms in half of the layers (e.g., 
Y-G layer) are replaced by Rh atoms with a Rh'^^ elec- 
tronic configuration, leaving 5 d-electrons in the outer 4d 
shell. Rh, being a lighter element, has a smaller atomic 
SOC than Ir. Hence, adding Rh atoms in this way cre- 
ates a staggered SOC along the c axis. We predict that 
tailoring a superlattice of Sr2lrRh06 will lift the nodal 
line present in SrIrOs and may produce a strong topolog- 
ical insulator in perovskite iridates. Sr2lrCo06 is another 
candidate as Co^+ in 3d shell leads to a similar effect of 
staggered SOC, but due to a smaller strength of the SOC 
in Co, stronger U may be necessary to generate a strong 
topological insulator. The band structure of the proposed 
materials and the effect of the combination of SOC and 
Hubbard interaction will be presented elsewhere.^ 

Material realization of topological insulators has so far 
been limited to narrow band gap semiconductors based 
on Hg or Bi, where s and p orbitals are involved. In 
addition to a recent theoretical proposal, where per- 
ovskite heterostructures of transition metal oxides grown 
along the (111) direction could be tuned to topological 
insulators^, our results further encourage the experimen- 
tal search for topological phases via synthesis and charac- 
terization of transition metal oxide materials. In particu- 
lar, the strong SOC in Ir makes iridates good candidates. 

In summary, we have shown that the underlying band 
structure of SrIrOa possesses a line node near the Fermi 
surface due to the crystal symmetry. Such a line node 
can be lifted by breaking the sublayer reflection symme- 
try and a transition to a strong topological insulator can 
be achieved. We suggest that superlattices of Sr2lrRh06 
or Sr2lrCo06 tailored in the direction suggested in this 
paper may be good candidates to achieve a strong topo- 
logical insulator. 
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